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Abstract 

Asymptotical behavior of the distribution function of local density of states 
(LDOS) in disordered metallic samples is studied with making use of the 
supersymmetric cr-model approach, in combination with the saddle-point 
method. The LDOS distribution is found to have the logarithmically nor- 
mal asymptotics for quasi-lD and 2D sample geometry. In the case of 
a quasi-lD sample, the result is confirmed by the exact solution. In 2D 
case a perfect agreement with an earlier renormalization group calcula- 
tion is found. In 3D the found asymptotics is of somewhat different type: 
P(p) ~ exp(— const | In 3 p\). 

PACS numbers: 71.55.Jv, 71.20.-b, 05.40.+j 
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I. INTRODUCTION 



Mesoscopic fluctuations of various physical quantities in disordered systems have been 
intensively investigated during the last years M . It was understood that a whole distribution 
function is to be studied to get the complete physical information concerning properties of 
a disordered system. In particular, statistical fluctuations of local quantities, such as eigen- 
function intensity and local density of states (LDOS) have attracted a considerable research 
interest recently. Distrbutions of eigenfunction amplitudes and of LDOS are relevant for de- 
scription of fluctuations of tunneling conductance across a quantum dot 0, of some atomic 
spectra properties ||, and of a shape of NMR line ||. The microwave cavity technique 
H allows to simulate a disordered electronic system and to observe experimentally spatial 
fluctuations of the wave intensity ||. 

Distribution of eigenfunction amplitudes \ipf{r)\, characterizes properly the spatial fluc- 
tuations in a closed system with well defined energy levels. On the other hand, in an open 
sample, connected to the leads, the states are broadened in an energy space. In this case, it 
is more appropriate to consider the distribution of LDOS, 

p(E,r) = -Tr-W^r, r; E) , 

where Gr is the retarded Green function. Distribution of LDOS in open metallic samples was 
studied by Altshuler, Kravtsov and Lerner |7]|| for the case of two- and (2 + e) -dimensional 
systems. The consideration was based on the derivation of the effective a-model and its 
subsequent renormalization group analysis, following the earlier ideas by Wegner |J. It 
was found that the LDOS distribution is close to the Gaussian one but has slowly decaying 
logarithmically normal (LN) asymptotics. When approaching the Anderson transition in 
(2 + e) dimensions, the distribution was found to cross over to the completely LN one. 



Altshuler and Prigodin [Kj studied the LDOS distribution in strictly ID chains and also 
found the LN form of the distribution. It was conjectured on the basis of this similarity 
|7|JT0[| that even in a metallic sample there is a finite probability to find "almost localized" 



eigenstates. 

More recently, it was recognized fll~l~|-|T4]j that the statistical properties of various quanti- 



ties characterizing fluctuations of the wave functions can be very efficiently studied with use 
of the supersymmetric approach. In particular, it was shown in [ |TT| , |Tlfl that the distribution 
P y {y) of eigenfunction amplitudes y = ^^(rol 2 (here V is the system volume introduced 
for normalization purposes) can be written in terms of the supersymmetric cr-model in the 
following way. Let us define the function Y(Q ) as 



Y(Qo)=[ DQ(r)exp{-S[Q}} , (1) 

where S[Q] is the cr-model action 

S {Q] = -Jd d rStr\^(\7Q) 2 -nvr ] AQ\ , (2) 



Q{r) is a 4 x 4 supermatrix field, A = diag{l, 1, —1, —1}, Str stands for the supertrace, 
D is the diffusion constant, rj is the level broadening (imaginary frequency) and v is the 



mean DOS fL9 |. We do not go into details of the supersymmetric formalism here, which can 
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be found e.g. in fT5HT7l. For the invariance reasons, the function Y(Qq) turns out to be 
dependent on the two scalar variables 1 < Ai < oo and — 1 < A2 < 1, which are eigenvalues 
of the "retarded- retarded" block of the matrix Q . Moreover, in the limit 77 — > (at a fixed 
finite value of the system volume) only the dependence on Ai persists: 



Y(Q )=Y{\ 1 ,\ 2 )^Y a (2nvri\ 1 ) 
The distribution of the eigenfunctions intensities is given by flIT,H 



(3) 



A 2 

dy' 
dy 2 



y (Ai 



y 



(4) 



The distribution function P p (p) of LDOS can be also expressed through the function 



y(Ai 3 A2) [ I3UIB |. For convenience, we normalize p by its mean value: p = p/(p), where 
(p) = v, and write again p instead of p henceforth. Then, 



P(p) 



Id 2 - 
- — ^—^ { 1 2 d\ 1 Y(X 1 




(5) 



where 



K(Ai 



1 AX 

CLA 2 -: 7— 

-1 Ai — A 2 



Let us note the important symmetry relation found in [jT3|: 



P(p- L ) = P 6 P{P) 



(6) 



(7) 



It follows from eq.(|5]) and is completely independent on a particular form of the function 
^(Ai,A2). Obviously, eq.(0) relates the small-p asymptotical behavior of the distribution 
P p (p) to its large-p asymptotics. 

Relations (f|) and (0) gave the possibility to study the distributions of eigenfunctions 
amplitudes and of LDOS in various regimes. Exact calculation of the eigenfunction ampli- 
tudes distribution in a quasi-lD sample of arbitrary length was performed in JTBj. LDOS 
distribution in a sample with given level width was calculated in [ p~2| , |T8f for the OD and 
pOj for the quasi-lD case. In Ref. [[L4|], the distributions were studied in the vicinity 
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of the Anderson metal-insulator transition, where Y(Q) acquires the meaning of an order 



parameter function. In Ref. |2l|j2^] the distribution P y (y) in a metallic sample was studied 
by a perturbative method which is valid for not too large y. 

In a recent paper, Muzykantskii and Khmelnitskii |23[ studied the long-time asymptotics 
of the current relaxation in a disordered metallic sample via the supersymmetric a-model 
approach combined with the saddle-point method. Fal'ko and Efetov [24| employed the 
saddle-point approach suggested in |2B| to study the asymptotical behavior of the distri- 
bution of eigenfunctions amplitude, P y (y), given by eq.(f|). In the present paper, we apply 
this method to investigate the asymptotical behavior of the distribution function of LDOS 
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in an open sample (i.e. with leads being attached). Precisely this distribution function was 
studied by Altshuler, Kravtsov and Lerner |7||| in the renormalization group approach, so 
that the direct comparison is possible. We find the LN asymptotics of P p (p) in 2 dimensions, 
in agreement with the result of @. As we discuss in Section 4, it is not completely obvious 
that this agreement should have been expected. 

In the case of a quasi-lD geometry, we also calculate the LDOS distribution function 
using the exact solution of the ID cx-model [E9. In the metallic regime, the LN asymptotics 
obtained by the saddle-point method is reproduced. In the deeply localized regime, the 
distribution function P p (p) takes a completely LN form. 



II. ASYMPTOTICAL BEHAVIOR OF THE LDOS DISTRIBUTION FUNCTION 
FROM THE SADDLE POINT APPROXIMATION. 

We consider a disordered sample in the weak localization regime, which means that 

fi' 1 d d q 1 

where £ is the usual parameter of the perturbation theory [0 . In this situation, the distribu- 
tion of normalized local DOS p is mostly cocncentrated in a narrow Gaussian peak around 
its mean value p — 1, with the width ((p — l) 2 ) ~ £ 1 This close-to-Gaussian shape 

of P p {p) holds however in the region \p— 1| <C 1 only. We will consider in contrast the "tails" 
of the distribution, p <C 1 and p 1, where a much slower decay of P p (p) will be found. 

As is seen from eq.(||), the LDOS distribution P p (p) is determined by the function 
Y(\i,\2) at Ai > Therefore, in the asymptotical region /) > 1, only the form of 

y(Ai,A2) at large values of Ai > p/2 ^> 1 is relevant. In this region, and under the con- 
dition (H), F(Ai, A2) can be found with exponential accuracy with use of the saddle-point 
method |23| , |24]| . We consider the case of an open samle and neglect any inelastic processes 
in the bulk of the sample. The energy levels are then broadened solely due to the possibility 
for a particle to escape to the leads, and the action is given by eq.(H) at 77 = 0, 

S[Q] = - />rStr^(VQ) 2 (9) 



The saddle point equations have the same form as in Ref. |23) for u = 0: 



A0i = A6 2 = 0, (10) 

where cosh#i = Ai; cos #2 — A2, and A is the Laplace operator. Assuming the perfect 

coupling to the leads, the corresponding boundary condidtion takes the form Q|] eac [ s = A, 
or in other words, 

0i, # 2 |i eads = (11) 

Another boundary condition which follows from the definition ([[]) of the function 

y(Ai°\A 2 °' ) ), fixes the value of 8\, 62 in the observation point ro, where the LDOS dis- 
tribution is studied: 
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cosh 0i (r ) = aS 0) = cosh^ 0) ; cos 6 2 (r ) = Af = cos^ 0) (12) 

We start from the case of a quasi-lD sample. Let the observation point be at x = 0, 
whereas the sample edges at x — — L_ and x = L + , so that the sample length is L = L_+L + . 
The solution of the saddle-point equation (fLOf) with the boundary conditions (0), (0) is 

i #i 0) fi - -T-) ; x > o 

^) = U 0) l w: _ as) 



^r(i-fr) ; x<o 

and similarly for 2 - The action on the saddle point solution is given by 

s = nr / l(v9l)2 + (V92)2] A = nr l ° r + * n (r + + z:) (14) 

We find therefore that in the region a?" » 1, 

F(A<»', Af) ~ e -» ~ exp {-^ (i- + J-) h,'(2Af ) j , (15) 

with exponential accuracy. Dependence of Y on Ajj is n °t important within this accuracy, 
since it gives simply a prefactor after the A2-integration in eq.@. Therefore, according to 
eq.(||) the asymptotics of the LDOS distribution function at p ^> 1 is 

P p {p) ~ F(A?» = p/2) ~ exp {-^ (i- + i_) ln 2 p} , (16) 

up to a preexponential factor. According to the symmetry relation ([?[), eq. (|T6|) is equally 
describes the small-p asymptotics of P p (p) at p <C 1. 

The found LN asymptotics of the LDOS distribution in an open sample should be con- 
trasted to the "square-root-exponential" behavior of P p (p) and P y (y) in a closed sample 
Tl| , p!7|p0| , |22| . Let us show how the difference appears in the present saddle-point approach. 



In a closed sample, the distribution of LDOS, P p (p) is defined only if a finite level width 77 
in the action, eq.(0), is introduced. It may account for some inelastic scattering processes 
in the bulk of the sample. The saddle point equation acquires then the same form, as in ref. 



DA9 1 - 2r7sinh0i = (17) 
Eq.(|l7D should be supplemented by the boundary conditions 

V ^illeads = ' ( 18 ) 



where V™ is the normal derivative, and by the condition ( |T2"D in the observation point. 

Let us first assume that the saddle-point solution satisfies the condition 9i(x) 3> 1 
throughout the sample. It will be seen below that this assumption is valid if E c /r) ^> 1, 
where E c = D/L 2 is the Thouless energy. Then sinh#i in eq.flTTj) can be approximated by 
e 6l /2, and the solution for x > is found to be 
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DC 



{l + cos [VC(L + -x)]} 



where the constant C is defined from the condition 



L+VC = 7r — arccos [ 1 



If A^ 0) > E c /ij, the solution (|T|), (GO) takes the form 



CD 



27/Ai' 



(0) 



(19) 



(20) 



The action 



D \ 71 



S ~ I dx 



1 — cos 



7T 



X + 



(0) 



(21) 



7T 



?i) 2 + vri/r/e 01 



(22) 



is then found to be dominated by the vicinity of the observation point, where eg. fl2l|) can 
be reduced to 



/ 



-0i(x) 



2X1 



(o) 



l + x> 



V 



(23) 



Substituting (p3~D in (|22|) and taking into account also an analogous contribution of the x < 
region, we get 



S ~ Birvy/DriXi 



,(°) 



In the opposite case 1 <C A^ <C i?cA7; the solution ([19|) , (|20|) is reduced to 



The corresponding action is 

5 ~ 2iiur]\ ( ] )) L 
Combining ([24]) and fl26|) , we find therefore for E c /i] ^> 1 

P P (p) ~ Y(^ } = P/2) ' 



exp{-47ri/V2£>77p} , p > ^ > 1 



exp{— TcvqLp} 



> p > 1 



(24) 



(25) 



(26) 



(27) 



If E c /r] < 1, one has to use the exact equation ( jT7[) , which has the following solution: 



Je 1 [coshtf-cosh^L)] 1 ^ 
where #i(£) is defined from the condidtion 



(28) 
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1/2 dd 



(29) 



i(L) [cosh-# — cosh6 l i(L)] 1 / 2 

We find that for E c < r\ the solution fl28|) in the vicinity of the observation point x <C \J D /rj 

has exactly the form ( p3|) for all 3> 1, so that the action is given by eq. (|24]) . Thus, in 
this case the first of the asymptotics ( p7|) is applicable 

E,. 



P p (p) ~ exp{-47rz/^2 J Dr/p} , — < 1 , p > 1 (30) 

Note that eq.(^) has exactly the same form as the asymptotics of P p (p) in the infifnitely long 
20|j . Eq.(|27j) is also in agreement with the asymptotical behavior of the distribution 
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function P y (y) of eigenfunctions intensities [ gj] : 



P y (y) ~ exp^^Wity/L} ; y » (31) 

which is immediately reproduced from eqs.(f|), (|3l"l) . 

We return now to the case of an open sample and consider 2D geometry. As was suggested 
in [p3~| , we suppose the sample to be a disk of a radius L surrounded by a perfectly conducting 
media, with the observation point located in the center of the disk. We will find the result to 
be logarithmically dependent on L, that justifies its application to a sample of more general 
2D geometry with a characteristic size L. 

Eq.([T0|) for 6 l 1 (r) with the boundary conditions (|TT|), ( |12| ) take now the form 

0''( r ) + = o (32) 
r 

#i(q = #i 0) (33) 

^(L)=0 (34) 

The boundary condition (^) has to be put not in the observation point r = 0, but rather 
at certain distance from it r = Z* This is because the cr-model approximation breaks 



down for momenta g > / 1 , where / is of order of the mean free path. We will specify the 
value of below. The solution of eq.(|3"2"|) with the boundary conditions (B3|), ([54]) reads 



0(0) 

0i (r) = 0? - ^j-jy In r/Z* , k < r < L , (35) 



or, for Ai(r) = cosh #i(r) ~ |e e ' l( ' r * ) 



1 ln-^L/i.) 

Ai(r) = AHr] (36) 




The condition of applicability of the cx-model (diffusive) approximation is [2^ d6i/dr < l/Z. 
It leads to the following equation for the cut-off scale Z* 

1 = Mi g{ 0) 1 , , 

I dr r=h \n{L/Q Z* ' 1 ' 
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The action (^]) is now given by 



S ~ /" d 2 r(V0i) 2 = 7i 2 vD rdrr[9[(r)f 



lnL/7* 

Therefore, the LDOS distribution has the asymptotics 



(38) 



f 7r 2 ZADln 2 p] , , 

P» ~ exp \~7TTT\ > ( 39 ) 



where Z* satisfies the condition 



ln(L/Z* 



= 'n^y (40) 



The present consideration is meaningful provided Z* <C L, which holds if In p <C L/l. Under 
this condidtion, eq.(f4Tj|) yields 

i* * ^-^rr ( 41 ) 



Taking into account the logarithmic dependence of the exponent in ( |3Tf ) on Z* and neglecting 
the corrections of the ln(ln. . .) type, we can approximate Z* in eq.(|3^) by I. The result is 
of the LN form, in exact agreement with the asymptotical behavior for P p (p) found by 
the renormalization group method in 0. We discuss this agreement in more detail in 
Conclusion. Let us also note that in contrast to the ID case, in 2D the asymptotical form 
of the distribution P p (p) in an open sample, eq. (p9[) , is similar to that of P y (y) in a closed 
sample |24j . 



Proceeding in the same way in 3D case, we get the equation 



0' (r) 

+ = > ( 42 ) 



with the same boundary conditions (|33|), (|34|). The solution has the form 

^(r)=^ 0) ^, U<r<L (43) 

The condition fl3~T|) for Z* now yields Z* ~ W^, and the action is estimated as 

S = 2ti 2 vD J drr 2 [e[(r)} 2 

= 1-K 2 vDUQf 2 
= 2tt 2 vDW[ 0)3 

= M¥i )3 , (44) 



8 



where g{l) = 2tt 2 uDI is the conductance of a cube of the size I. Therefore, the LDOS 
distribution asymptotics is 

P p (p)~exp{-s(/)ln 3 p}, (45) 

Let us note that I is of order of mean free path, so that it is defined up to a numerical 
coefficient of order unity. This ambiguity does not affect the leading order behavior of 
the exponent in eq.(|3^) for the 2D case, due to its logarythmic dependence on /. On the 
other hand, in 3D case, the exponent in ([EJ) is proportional to I. To fix the corresponding 
numerical coefficient, one should go beyond the long- wave- length a-model approximation 
and consider the problem in the region r < 

To conclude this section, we have shown that the LDOS distribution P p (p) has logarith- 
mically normal asymptotics (|i~6D , (j39|) at p ^> 1 and p C 1 for quasi-lD and 2D samples, 
whereas for 3D systems a somewhat different result (^) was found. In the next section, we 
confirm this result for the quasi- ID case by the exact solution of the problem. 



III. LDOS DISTRIBUTION IN A QUASI ID SAMPLE: EXACT SOLUTION. 

In the case of a quasi-lD geometry of a sample, evaluation of the functional integral 
in ([]]) can be reduced to the solution of an evolution equation of the diffusion type for the 
function Yi( A i, A 2 ; x) |i~5| : 



ox 



I A T— 1 - - 

J del de\ + lW 2 W 2 



Tfo a \ |sinh 6 l sin 9 2 \ ( 
(cosht^ — cos 62) 

For the general case of a finite quasi-lD sample with distances L_, L + from the observation 
point to the edges, the function y(Ai,A2) defined in the preceding section is given by the 
product 

y(Ai, A 2 ) = Yi(Ai, A 2 ; L_)F 1 (A 1 , A 2 ; L+) (47) 

For simplicity, we will concentrate on the case of a semi-infinite sample L + = 00, L_ = L. 
In this case F(Ai, A 2 ) = Y\(\i, A 2 ; L), and we will omit the subscript "1" in the notation 
"Yi". Results for the general case will be presented in the end of the section. Solution 
of eq. fl46|) with the boundary condition corresponding to a perfectly conducting lead was 
recently found to be 

00 roc 

Y{Xx, A 2 ; L) = 1 - (A x - A 2 ) / dkc mk P m {X 2 )P- 1/2+ik/2 {X 1 )e- e ^ t ; 

t= 2^D- ^ - (2m (tl )k + ipTF^ ; ^ = (»+V2) 2 + *V4 (48) 
Substituting eq. (|48|) in (|), we find 
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Y(X 1 ) = In ^l±i - 2 jH ^^A; tanh(vrfc/2)P_ 1/ 2+ife/2 (A 1 ) e -'( fe2+1 )/ 4 (49) 

The first term in eqs. (fi8|) and ( |49|) corresponds to the limit of a closed system with no level 
broadening, in which case Y(Xi, A 2 ) = 1 and P p (p) = S(p). It can be really proven by direct 
substitution of the first term of eq. (f49|) into eq.([5]) that the resulting contribution to P p (p) is 
equal to zero at any p ^ 0. To evaluate the second term, we use the asymptotic expression 



for the Legendre functions P_x/ 2 +i 7 (£) at z ^> 1 |26 



P - 1/2+ ^ } -v^l r(i/2 + . 7 ) +e r(i/2-fy)/ (50) 

This allows us to reduce the second term in (f49[) to the form 



At In Ai ^> 1 the integral can be evaluated via the saddle point method. The saddle point 
is k = i ln2Xl and yields the following contribution: 



2 1 


r 


i 


1 ln2Ai\ 
' 2t J 


tAi i /"ln2Ai^ 


2 r 


f In2Ai ^ 

v 2t ; 





i ln 2 2Ai 



Fs.pXA,) * -4W— - — , ^ exp -- - — (52) 



4 4t 

In fact, when deriving eq. Q52fl , we should shift the integration contour in the complex plane 
in order that it would pass through the saddle point. However, if ln2Ai > t, the contour 
crosses then the pole of the integrand, k = i. Evaluating this pole contribution, we find that 
in cancels exactly the first term in eq.flHf), (|49"]). Thus, we get 



F (Ax) ~ In (In 2Ai -t) + F s .p. (Ax) (53) 

where -&(x) is the step function and F s .p. is given by eq. (|52|) . Now we substitute this result 
into the formula (HD for the LDOS distribution. The leading contribution at p ^> 1 is found 
to be given by the second term, Fs.p. (Ax), in (|53|). It can be found by noticing that the 
factor 





1 


r 


3 


1 ln2Aj\ 
+ 2i j 


1 - 




2 r 


( ln2Ai \ 

V 2t ; 





varies logarithmically slow with Ax, and can be simply taken at Ax = p/2. The remaining 
integral is of the form 



1 ft ln 2 2Ax 
a Ax exp <^ 

p/2 y/X^-p/2) I 4 



it 



dy it In p In p In y In y j 

y^T) eXP 1^4 4t 2t 4TJ 

B 2-if H -4— ir (54) 



10 



where B(x,y) is the Euler's beta function. When writing the last line, we assumed that 
In p > y/i. It is not a restriction in the metallic (short sample) case, t < 1, and includes the 
region of characteristic LDOS, hip ~ t, in the insulating regime, t 3> 1. 

Collecting now all the remaining factors in eqs.(||), (|52D, we finally get the distribution 
of LDOS: 

P p (p) ~ -^ eX p(-i(t + lnp) 2 l (55) 



2pV7rf I it 

Let us remind that eq.(^) has been derived under the assumption p ^> 1. An additional 
condition In p ^> \ft in the insulating regime is not very restrictive, since the characteristic 
scale of In p in this case is In p ~ t, as is seen from (f^). Moreover, calculating the integral 
([54]) in the opposite case lnp -C t, we arrive at the same LN behavior as in eq. (|55|) ; only the 
pref actor is different: 

W^ ! ^exp{-l(t + lnp) 2 } (56) 

Therefore, we have found the LN asymptotic behavior (p5|), ( |56"D of the LDOS distribution 
P p (p) ar /3 > 1. The symmetry relation fl7|) allows us to extend its validity to the region 
of p 1 as well. Note that eqs.(|55|, ([Blf) completely preserve their form under the trans- 
formation (^). Eqs.(^), (^) are in full agreement with the results of the saddle point 
approximation of section 2. 

Let us remind theat the above results have been derived for the case of a semiinfinite 
sample: t + = oo, t_ = t, where t± = L±/2irvD . Now we consider briefly the general 
case of finite t + , t_. In this case the function K(Ai, A2) determining the LDOS distribution 
is defined by eq.(|47|) with Y"i(Ai, A2; L) given by eq.(^8|). Using the above analysis of the 
asymptotic behavior of Yi at Ai ^> 1, we find with an additional condition ln2Ai > t + ,t_, 

F(A i; t+,t_) = yF s .p.(Ai,t + )F s .p.(Ai,t-) , (57) 

with V s .p.(Ai,t) given by eq. Q52]) . This leads to the following result for the LDOS distribu- 
tion: 



In p 
2t7 




P p (p) ~ . — - exp + 1_) - - - + - \\n 2 p (58) 




where 



(x - 1/2) 



(x + l/2)r(ar) 

Eq. (|58D holds in the metallic regime t+,t- < 1 for all p ^> 1 and is again in agreement with 
the result of the previous section, eq . (|T6|) . In the insulating regime eq. (|58|) holds provided the 
additional condition lnp > is satisfied and represents therefore the far asymptotics. 

In this case, there is also an intermediate regime 1 < lnp < t + , where we assumed for 
definitness t- < t + . Then K(Ai, ^2,t+) can be approximated by 1, and we find P p (p) given 
by eq.(|55|) with t = i_. In all the cases, the equality (|^) allows to get the asymptotical 



behavior in the region of small LDOS, p <C 1. 
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In conclusion of this section, let us note that in the insulating regime, t-,t+ 3> 1, the 
obtained results for P p (p) are completely analogous to what was found by Altshuler and 
Prigodin [|10j for the case of a strictly ID sample by Berezinskii technique. This confirms 
once more the general conjecture Jl7j that the statistics of smooth envelopes of the wave 
functions in ID and quasi-lD samples are equivalent. 



IV. DISCUSSION AND CONCLUSION. 

In this article, we have studied the asymptotical behavior of the LDOS distribution 
function P p (p) in disordered metallic samples at p <C 1 and p ^> 1. For this purpose, we 
used the supermatrix a-model approach, which allows to express the distribution P p (p) in 
terms of the function Y(X 1 , X 2 ) defined as a certain integral over the supermatrix field. In 
the asymptotical regime, p C 1 and p > 1, this integral can be estimated via the saddle- 
point method, leading to the log-normal asymptotics of P p (p) for quasi-lD and 2D sample 



geometries, and to the result (|5|) in the 3D case. 

In quasi-lD the result can be also found from the exact solution, which indeed yields 
the LN asymptotical behavior in the case of a short (metallic) sample. In the opposite case 
of a sample much longer than the localization length, the whole distribution function P p (p) 
has the log-normal form. This is completely analogous to the form of the LDOS distribution 
in a strictly ID sample in the strongly localized regime studied in |10| by the Berezinskii 
technique. 

In 2D, the obtained asymptotics of the LDOS distribution is in full agreement with the 
result of renormalization group treatment M. This agreement is highly non-trivial, for the 
following reason. The RG treatment is based on a resummation of the perturbation theory 
expansion and can be equally well performed within the replica (bosonic or fermionic) or 
supersymmetric formalism. At the same time, the present approach based on the super- 
symmetric formalism relies heavily on the topology of the saddle-point manifold combining 
non-compact (Ai) and compact (A 2 ) degrees of freedom. The asymptotic behavior of P p (p) 
at p ^> 1 and p <C 1 is determined by the region Ai ^> 1 which is very far from the "per- 
turbative" region of the manifold Q ~ A (i.e. Ai, A2 — 1). It is well known [E7| that for the 
problem of energy level correlation, the replica approach is not able to reproduce the correct 
result (which can be obtained by the supersymmetry method) even if all orders of perturba- 
tion theory are taken into account. This is related to the fact that within the replica method 
the topology of the saddle-point manifold is not properly reflected. This might suggest a 
conclusion that any result based on a resummation of the perturbative expansion (which 
can be therefore obtained within the replica approach) is incorrect for the same reason. 
The agreement of the results of supersymmetric and renormaliztion group treatments of the 
LDOS distribution shows however that this conclusion would be wrong. Apparently, this 
problem is of a diferent kind, as compared to the level correlation one, so that the replica 
trick combined with the RG is able to reproduce the correct result. The difference seems 
to be the following: in the case of level correlator both compact and non-compact sectors 
in the supersymetric formulation are equally important, whereas for the present problem 
only the non-compact sector was essential, with compact one playing an auxiliary role. Let 
us note that the same situation appears in the vicinity of the Anderson transition f28|JH}| 
where the function y(Ai, A2) acquires a role of the order parameter function and dependes 
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on the non-compact variable Ai only. The above agreement obtained for the LDOS dis- 
tribution provides therefore support to other results obtained with making use of the 
renormalization group approach and 2 + e expansion. 
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